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ABSTRACT 

The  results  and  findings  of  researches  into  stress  wave  propagation 
of  granular  media  are  presented.  The  purpose  of  this  work  is  to  clarify 
air- induced  ground  shock  phenomena  in  soils,  leading  to  reliable  analyti¬ 
cal  methods  of  prediction,  and  to  provide  guidance  for  experimental  veri¬ 
fication  and  research. 

This  final  report  summarizes  work  on  one- dimensional  wave 
propagation  of  bi- linear,  locking,  and  compacting  materials,  and  presents 
the  first  results  of  the  extension  of  these  researches  into  two-dimensional 
wave  phenomena  in  such  non-linear  media. 
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I.  Introduction 


This  report  Is  a  G\itmnax*y  outline  of  the  work  vhlch  has  been  ac- 
complished  under  our  current  contract,  AF29(6oi)-2855»  the  purpose  of 
Tihich  vas  to  investigate  stress  vave  propagation  phenomena  in  granular 
materials.  Detailed  results  of  this  vork  have  been  issued  in  three 
technical  reports; 

"Attenuation  of  Stress  Waves  in  Bi-Linear  Materials", 
by  Richard  Slcalak  and  Paxil  Weidlinger,  AFSWC-TN-60-30j 
October  i960. 

"A  Method  for  Prediction  of  Groxind  Shock  Phenomena 
in  Soils"  (Secret)  by  Paul  Weidlinger  and  Alva  Matthews, 
ArswC-TDR-6l-66i  August  I96I. 

"step  Load  Moving  with  Low  Subseismic  Velocity  on  the 
Surface  of  a  Half -Space  of  Granular  Material",  by 
Hans  H.  Blelch  and  Ewald  Heer,  AFSWC-TDR-63-2;  April 

1963. 

This  work  is  a  continuation  of  a  broad  effort  on  this  subject,  and 
previovis  results  have  been  published  (see  References  1,  5,  o). 

This  research  effort  in  its  narrower  sense  is  directed  towards 
the  clarification  of  air-induced  groxond  shock,  but  it  also  has  a  direct 
bearing  on  some  other  cxrrrent  investigations  relating  to  the  effect  of 
direct  groxind  shock. 

Stress  vave  phenomena  in  granular  materials  open  up  an  entirely 
new  field  of  Investigation  where  scarcity  of  both  exi>erlmental  and  ana¬ 
lytical  resxxlts  is  prevalent.  The  earliest  investigations  by  our  group 
date  back  to  1958  (Reference  ll)  and  to  195^  by  USSR  researchers  (Refer¬ 
ence  15). 
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The  ultimate  purpose  of  this  work  is  to  malte  available  methods 
which  will  be  suitable  for  the  design  of  hardened  underground  construc¬ 
tion;  the  specific  purpose  is  to  provide  sufficient  theoretical  back¬ 
ground  for  the  performance  and  guidance  of  laboratory  and  field  experi¬ 
ments. 

In  investigations  which  are  restricted  to  one  dimensional  pleuie 
wave  phenomena,  a  simple  stress -strain  relation  will  describe  the  be¬ 
havior  of  the  medium  if  strain  rate  effects  can  be  neglected.  Experi¬ 
mental  evidence  indicates  that  under  one  dimensional  dynamic  compression, 
gramHar  materieils  show  a  hardening  stress -strain  relation^  that  is,  the 
second  derivative  of  the  stress-strain  cvirve  is  positive.  This  type  of 
stress-strain  curve  has  been  approximated  by  various  slrngpllfied  models, 
such  as  locking  materials  aixd  bi-linear  materials.  It  wgis  also  shown 
that  closed  form  or  numerical  solutions  can  be  obtaj.ned  for  one  dimen¬ 
sional  waves  in  compacting  medium,  which  has  a  more  gexieralized  form  of 
stress-strain  relation,  and  which  is  a  good  approximation  of  some  experi- 
mentedly  obtained  stress-strain  curves. 

As  we  proceed  into  investigation  of  cases  where  more  than  one 
spatial  dimension  is  involved  or  where  the  wave  front  cannot  be  con¬ 
sidered  plane,  information  in  addition  to  the  stress-strain  curve  is  re¬ 
quired.  The  case  of  comi)lete  spherical  symmetry  which  can  be  treated 
mathematically  as  a  one  dimensional  case  requires,  for  Instance,  the  in¬ 
troduction  of  a  yield  condition  which  in  granular  materials  can  be  de¬ 
rived  frcxn  the  assunption  that  Coulomb's  law  of  failure  is  applicable 
under  dynamic  conditions.  A  separate  law  assuming  a  locking  property  of 
the  medium  must  be  postulated  at  the  wave  front. 


As  ve  Investigate  general  tvo  dimensional  caaes  in  addition -to 
the  yield  condition,  other  properties  of  the  medium  must  be  postulated^ 
such  as  flotf  rule,  which  can  be  expressed  in  relatively  simple  terms  in 
the  cctse  of  an  isotropic  material.  It  turns  out  that  in  addition  to 
the  flow  rule,  at  least  one  further  relation  is  required,  which  can  be 
obtained  by  further  postulating  the  existence  of  a  plastic  potential. 

means  of  these  assumptions,  the  specific  problem  of  stress 
waves  generated  by  a  step  load  moving  with  low  subselsaic  velocity  was 
ccxtsidered  and  solved  in  our  third  report  under  the  current  contract. 

In  the  investigation  of  this  problem,  numerous  assumptions  regarding 
the  behavior  of  the  medium  can  be  made  which  yield  solutions  which  ai>- 
pear  to  be  physically  significant.  Each  of  these  solutions,  which  spe> 
clflcally  consider  whether  the  medium  is  con^esslble  or  IncooqareBsibla, 
resulf.f  in  a  slightly  different  behavior  of  the  material.  Further  pet- 
slble  modes  of  response  can  be  postulated  as  other  cases  and  other  velo¬ 
cities  of  the  moving  step  load  are  considered.  The  idausibllity  of  some 
of  these  ass^ill^tlons  can  be  suppoirted  by  the  conversion  of  the  mathe¬ 
matical  model  into  a  sliqpllfied  physical  mechanism  which  in  some  re¬ 
spects  can  be  considered  analogous  to  the  sprlng-dashpot  model  of  stan¬ 
dard  solids.  Nevertheless,  the  final  choice  between  the  various  mathe¬ 
matically  admissible  modes  of  behavior  must  ultimately  be  resolved  by 
laboratory  and  field  e3q>erimeat8.  Consequently,  the  analytical  work 
Itself  is  directed  towards  the  solution  of  specific  problems  which  are 
reproducible  under  laboratory  conditions  or  in  field  experiments.  The 
results  so  far  obtained  have  shown  that  some  of  the  quantitative  results 
(such  as  attenuation  of  stresses  and  pairtlcle  velocities)  show  close 


agreement  vlth  field  measvirementa  of  weapons  tests^  and  the  existence 
of  a  dust  precursor^  (which  was  observed  in  weapons  tests)  hats  been 
demmstrated.  These  facts  are  liqportant  to  the  extent  that  they  Indi¬ 
cate  that  the  basic  assvoQitlonB  are  plausible  and  are  suppoirted  by  ex- 
Xjerlmental  evidence.  On  the  other  hand,  because  of  the  great  difficul¬ 
ties  which  are  inherent  in  laboratory  work  on  granular  materials,  other 
significant  phenoeana  can  oxily  be  discovered  (or  distinguished  fron 
secondary  effects  and  noise)  with  the  aid  and  guidance  of  quantitative 
predictions  based  on  the  previously  dlsctissed  theoretical  analysis. 

It  is  hoped  that  the  current  work  will,  in  fact,  lead  to  certain 
definite  statements  and  clarification  regarding  the  behavior  of  real 
materials  through  the  selection  of  the  pax^icular  mathematical  model  or 
mode  of  response,  which  corresponds  to  the  results  of  experimental 
measurements.  It  is  further  anticipated  that  tvo-dlmexislonal  investi¬ 
gations  will  make  it  possible  to  ccmduct  investigations  into  diffrac¬ 
tion  phezKxaena  which  are  1  implicit  in  laboratory  experiments  due  to  the 
effect  of  ffleasxuring  devices  and  the  finiteness  of  the  space  in  which 
the  ejqperlment  is  performed.  Finally,  the  solution  of  the  diffraction 
problem  can  leed  to  the  determination  of  the  effects  of  interaction  of 
burled  structures  with  stress  waves  in  granular  soils. 


II.  One  Dimensional  Plane  Waves 


One  dimensional  plane  stress  waves  are  characterized  hy  a  plane 
wave  front  with  zero  strain  gradient  parallel  to  the  front  and  generally 
a  non-zero  gradient  perpendlctilar  to  the  front. 

The  medium  ahead  of  the  front  is  at  rest  and  the  stress  is  e<]ual 
to  zero.  Behind  the  front  the  medium  is  in  a  state  of  plain  strain. 

Fbr  this  reason;  the  relevant  physical  characteristics  of  the  medium 
are  sufficiently  defined  hy  the  initial  density  and  the  dynamic  stress- 
strain  law  under  uniaxial  confined  compression.  Ilie  simplicity  of  the 
wave  phenomena  xmder  these  assumptions  admf.p  a  variety  of  simple  solu¬ 
tions  which  depend  only  on  the  choice  of  a  few  parameters  defining  the 
physical  characteristics  of  the  medium  and  that  of  the  applied  pressure 
pulse. 

Such  one  dimensional  plane  waves  are  generated,  if  a  saai-infinite 
medium  is  loaded  by  a  pressxire  pulse  of  \xnlform  intensity  over  the  en¬ 
tire  free  boundary.  One  dimensional  analysis  will  also  be  approximately 
valid  if  the  free  boundary  is  loaded  by  a  pulse  which  advances  over  the 
surface  with  a  velocity  which  is  substantially  greater  than  the  propa^ 
gation  velocity  in  the  medixim.  Hiis  type  of  loading  is  characteristic 
of  the  blast  pressure  generated  by  a  nuclear  ejqolosion  and  consequently 
resxilts  derived  from  plane  wave  analysis  can  be  used  to  provide  infor- 
znation  regarding  the  effects  of  air-induced  ground  shock  in  the  supei'- 
seismic  range,  which  in  many  soils  corresponds  to  peeik  asijient  pressure 
above  the  100  psi  contours. 


The  simplified  analysis  also  assuxoes  that  the  soil  Is  honogeneous 
and  temperature  effects,  phase  change^  or  transitions  can  be  neglected. 
Subject  to  these  limitations  the  results  of  analytical  predictions  are 
experimentally  verifiable,  provided  that  other  relevant  characteristics 
of  the  medium,  i.e.  stress-strain  law  auid  initial  density  are  known. 

1.  Stress -Strain  Characteristics  of  Granxilar  Soil 

The  classical  literature  on  soil  mechanics  provides  extensive  in¬ 
formation  regarding  the  static  stress-strain  characteristics  of  soils. 
More  recently  such  ejqperitnental  work  has  been  extended  to  the  dynamic 
behavior  of  fully  confined  soil  samples  (References  1,  2,  3,  9).  These 
ejqperiments  generally  seem  to  agree  on  the  following  points: 

a)  The  material  shows  an  initial  elastic  and  possibly  vlsco-elaa- 
tic  behavior  at  low  stress  levels  and  low  density. 

b)  At  higher  stress  levels  large  irrecoverable  strains  are  mani¬ 
fested  together  with  pleetic  behavior. 

c)  With  increeislng  pressvures  a  marked  decrease  of  compressibility, 
i.e.  a  hardening  in  the  stress-strain  curves  is  observed  and  unloading 
occurs  on  a  nearly  vertical  breuich  with  large  irrecoverable  strains 
accompanied  by  a  significant  increase  in  density. 

d)  Upon  reloading; the  material  follows  the  unloading  curve  up  to 
the  original  stress  level  and  from  there  on  shifts  into  the  original 
loading  branch  of  the  diagram.  At  these  hi^  stress  levels  no  signi¬ 
ficant  strain  rate  effects  are  observed. 

The  typical  form  of  such  a  stress-strain  diagram  is  shown  in  Fig¬ 
ure  1.  It  also  should  be  noted  that  in  some  materials  the  initial 


elastic -pleistic  behavior  is  not  observed  and  the  effect  of  reduced 
conqireBsibility  (i.e.  hardening)  is  the  governing  phenomenon.  At  the 
present  state  of  the  art  no  direct  relationship  between  the  shape  of 
the  stress-strain  curve  and  other  physical  properties  of  the  material, 
such  as  grain  size,  grain  size  distribution,  cohesion,  etc.,  is  estab¬ 
lished. 

2.  Idealized  Stress-Strain  Belations 

To  aJLlov  an  analytical  treatment  of  one -dimensional  soil  dynamics 
problems^ various  stress-strain  relations  may  be  assumed.  The  validity 
and  applicability  of  these  simplified  stress-strain  laws  depend  on 
pressvire  levels  and  on  the  purpose  for  which  these  analytical  models 
are  to  be  used. 

In  the  past  several  investigations  have  been  conducted  which  are 
based  on  linear  elastic  behavior.  Visco  elastic  behavior  of  a  modi¬ 
fied  Voigt  solid,  represented  by  a  spring  dashpot  model,  was  also  con¬ 
sidered  (Reference  4). 

The  research  which  is  reported  under  this  contract  is  based  on 
veurious  non-linear  materials  of  limited  compressibility,  the  sinplest 
of  which  is  a  stress -strain  relation  known  as  the  ideal  locking  medium. 
(See  Figure  2).  The  stress-strain  law  of  the  material  is  such  that  it 
exhibits  no  resistance  up  to  a  critical  strain,  which  corresponds  to 
the  theoretical  limit  of  its  compressibility.  Since  the  medium  cannot 
be  conqpressed  beyond  this  limit  the  stress-strain  diagram  becomes  verti- 
csLl.i.e.  paredlel  to  the  stress  axis  both  on  loading  and  unloading.  Other 
variations  of  this  locking  medium  have  eLLso  been  considered  which 


assume  a  different  type  of  response  before  the  critical  strain  is  .reached^ 
such  as  an  elastic  behavior  (Figure  3)  or  an  eleistic-plastic  behavior 
(Figure  4).  Stress  wave  phenomena  in  materials  of  these  types  have  been 
treated  in  Reference  5«  Spherical  wave  propagation  heis  also  been  con¬ 
sidered  in  Reference  6. 

Another  group  of  models  proposes  a  bi-lineal  stress-strain  re¬ 
lation.  In  this  medium,  in  its  simplest  form,  no  resistance  up  to  a 
critical  strain  is  assumed  euic,  beyond  this  point,  linear  elastic  re¬ 
sponse  is  postulated  (Figure  5).  A  different  assumption  regarding  the 
medium  before  it  reaches  the  critical  strain  leads  to  a  bi-linear  mater¬ 
ial  which  shows  an  initial  but  smaller  elastic  resistance  before  the 
critical  strain  is  reached  (Figure  6).  The  proi)erties  of  such  media 
were  analyzed  in  our  first  reiwrt  TII60-30  under  our  cunrent  contract 
(Reference  T).  By  varying  the  elastic  modulus  of  the  second  branch  of 
the  bi-linear  model  the  entire  domain  of  the  linear  elastic  and  the 
elastic -locking  media  is  included  (Figure  T).  Wave  propagation  pheno¬ 
mena  in  the  bi-linear  model  have  been  treated  in  Reference  T  for  a  dis¬ 
continuous  pressure  input,  i.e.  a  suddenly  applied  load.  The  formation 
and  propagation  of  the  strong  shock  for  a  continuous  input  hav'^  ‘  t. 
treated  recently  in  Reference  8. 

A  special  case  of  bi-linear  medium  is  also  that  of  a  material 
which  exhibits  linear  stress-strain  relation  upon  loading  and  follows 
a  steeper,  but  also  linear  branch  upon  unloading.  The  unloading  branch 
may  become  completely  vertical  approaching  the  behavior  upon  unloading 
of  that  of  a  locking  material  (see  Reference  5  end  Figure  8). 


All  the  above-deocribed  etresB-Btrain  diagrama  repreeent  an  at- 
ten5)t  to  approximate  the  actual  behavior  of  granular  soils,  i.e.  they 
represent  an  idealized  model  of  the  actual  stress-strain  diagrams  shown 
in  Figure  1. 

Numerical  calculations  have  shown  that  beyond  a  definable  depth 

from  the  surface  the  decay  of  the  peak  stresses  and  particle  velocities 
is  not  too  sensitive  to  the  general  shape  of  the  stress-strain  curve 
as  long  as  a  general  liardening,  that  is,  limited  compressibility,  is 
assumed.  In  investigations  which  are  directed  toward  free  field  phenomena 
at  a  great  distance  from  the  surface^  the  distinction  between  various 
types  of  bi -linear  and/or  locking  behavior  is  only  of  limited  importance 
(Reference  7)* 

In  the  second  report,  TDn6l-66  under  our  current  contract  (Re¬ 
ference  9i  an  attempt  has  been  mode  to  consider  e  more  general  form  of 
strain  hardening.  In  this  material  the  stress-strain  curve  bends  to¬ 
wards  the  a  axis  with  increasing  intensity  and  upon  uiHoadlng  shows 
Irrecoverable  strains  at  each  stress  levels  that  Is,  the  unloading 
branch  is  assumed  to  be  vertical.  Specifically,  It  Is  assumed  that 
the  stress-strain  curve  Is  Independent  of  strain  rate  effects.  Upon 
loading  It  is  characterized  by 


0 


and  upon  unloading  it  follows  the  vesrtlcal  line  defined  by 


t 


e 


max 
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It  reloads  on  the  unloading  curve  up  to  the  maximum  stress  and  beyond 
that  it  follows  the  original  loading  curve.  Since  the  medium  is  later¬ 
ally  confined,  a  constant  and  determinable  lateral  confining  pressure 
exists,  while  zero  lateral  strain  implies  that  the  vertical  strain  is 
defined  by 


e 


where  is  the  initial  density  (at  o  =  O)  and  the  density  increases  to 
p  at  the  strain  e.  It  is  further  assxuaed  that  the  material  has  no  ten¬ 
sile  strength.  A  medium  which  has  these  characteristics  is  called  a 
compacting  mediiom  (Figure  9)*  The  entire  domain  of  such  materials  is 
boiinded  by  two  extreme  forms  of  the  stress-strain  law.  In  one  extreme 
we  have  the  linear  material  shown  in  Figure  8  and  in  the  other  extreme 
the  ideal  locking  medium  shown  on  Figure  2. 


3.  One  Dimensional  Stress  Waves  in  a  Bi-I<inear  Material 

Consider  a  semi -infinite  speuie  characterized  by  a  stress -strain 
curve  shown  on  Figure  5  which  is  sxibjected  to  a  discontinuous  pressure 
pulse  over  its  free  boundary.  The  shock  velocity  z  will  be 


where  E  is  the  secant  modulus  corresponding  to  the  stress  level  a  at 
s 

the  front.  The  one  dimensional  wave  equations  in  Lagrangian  coordinates 


d^o  2 

i?  '  ' 


are 


2  o^u 
°  dz^ 

vhere  u(Z|t)  is  the  displacement  of  a  particle  at  a  time  t  from  its  Ini- 

2  Eti 

tlal  i)ositlon  z,  c  =  ^  ,  and  is  the  tangent  modulus  of  elasticity. 

These  vave  equations  can  be  solved  by  the  numerical  technique 
outlined  in  Reference  7  and  they  permit  the  evaluation  of  the  influence 
of  the  parameters  of  the  medivun  on  the  free  field  quantities.  One  of 
the  significant  findings  is  that  beyond  a  predictable  depth  the  peak 
stress  and  peak  particle  velocity  of  a  bi-linear  medium  become  nearly 
coincidental  with  that  of  an  ideal  locking  medium  of  Identical  critical 
strain  .  This  particular  depth  corresponds  to  the  time  at  which  the 
first  and  principal  rarefaction  wave  from  the  surface  arrives  at  the 
shock  front.  The  practical  implications  of  this  are  that  the  attenua¬ 
tion  of  stress  waves  in  the  far  field  depends  principally  upon  the  value 
of  the  critical  strain  and  the  decay  time  of  the  applied  pulse. 


One  Dimensional  Stress  Waves  in  a  Compacting  tfedi^m 

Consider  a  semi -infinite  space  characterized  by  a  stress -strain 
curve  of  the  type  shown  in  Figure  9»  subjected  to  a  discontinuous  pres¬ 
sure  pulse  p(t).  The  shock  velocity  of  z  will  be 


and  the  equation  of  motion 


is  defined  by 

dOj^/dt 

— - -  \  p  dt  -  P  - 

Cf„  J  z 

Z  o 


where  is  the  stress  at  the  front. 
2 
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We  note  that  in  the  case  of  =  E,j,  the  above  equation  degenerates  to 
that  of  the  linear  material  of  Figure  8  and  in  the  case  of  i;j,  =  “  to 
that  of  the  ideal  locking  medium  of  Figure  2.  The  general  form  of  the 
equation  can  be  solved  by  numerical  techniques,  but  a  useful  closed  form 
solution  is  obtainable  for  the  case 


=  constant 


which  isiplies  a  parabolic  stress-strain  curve  of  the  form 


0  .  (|)' 


where  k  >  0  and  n  ^  1  are  constant  characteristics  of  the  medium.  The 
case  n  =  1  implies  a  linear  medium  and  n  =  •  implies  the  locking  medium. 
Particular  values  of  n  and  k  can  be  selected  to  approxLinate  experimen¬ 
tally  determined  stress-strain  curves.  The  intensity  of  the  stress  at 
the  front  is 


Tt  n-l' 


where  I  =  j  p  dt  and  the  front  location  is  given  by 


n+1 
/  2n  vTrT 


t  n-l  i 

i 


The  above  equationa  define  the  attenuation  of  the  peak  stress  with  depth 
and  permit  the  evaluation  of  all  other  free  field  quantities. 

5.  Application  to  the  Prediction  of  Ground  Shock  Effects 

Of  the  various  one  dimensional  theories^  that  of  the  compactixig 
medium  turns  out  to  be  particularly  suitable  for  the  prediction  of  stress 
wave  phenomena  in  soils.  Currently  under  various  contracts  laboratory 

and  BE  field  experiments  are  in  progress  in  testing  the  applicability  of 
this  theory. 

In  our  second  report,  TDR6I-66  (Reference  9]i  it  vaa  shown  that 
ground  motion  data  obtained  for  nuclear  field  tests  can  be  reproduced 
by  the  application  of  this  theory  and  it  was  also  shown  that  the  i>re> 
diction  of  alr«lnduced  ground  shock  effects  cozresponds  fairly  closely 
to  most  currently  proposed  es^lrlcal  prediction  procedures.  The  corn* 
pactlng  medium  theory  seems  also  successful  in  defining  shock  spectra 
within  frequency  raises  of  ciurent  practical  interest. 

In  Table  1  below  a  comparison  of  quantities  obtained  by  use  of 
the  compacting  medium  theory  is  shown  with  the  semi -empirical  methods 
proposed  in  Reference  10. 
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TABLE  I 


Comparison  of  Theory  vith  ASCE  Manual  Predictions 


Weapon:  W  =  8  JIT 

Soil:  V  =  115  #/ft3 


P  =  200  psl 
o 

c  =  2000  ft/sec 


Maxlmuii  Effects  at  Depth,  in  Feet 
Theory  of  this  Paper  ASCE  Manual  Predictione 


Quantity 

O' 

50' 

100' 

O' 

50' 

100' 

Peak  Stress 

200  psi 

l84  psi 

169  psi 

200  psi 

180  psi 

160  psi 

Peak  Particle 
Velocity 

Vertical 

4.36 

sec 

4.01  — 
sec 

3.69-^ 

sec 

4.16 

sec 

3.75  — 

'  sec 

3.33  i!L 
sec 

Horizontal 

2.52  '■ 

2.32  " 

2.13  " 

2. 75  " 

2.5 

2.16  ’’ 

Displacement 

Vertical 

12.3  in 

10.9  in 

9.6  in 

14.5  in 

13.3  in 

12.0  in 

Horizontal 

Y.l  " 

6.3  " 

5.5  " 

00 

4.4  " 

4.0  " 

Ill .  Two  Dimensional  Problems  of  Wave  Propagation 


The  knowledge  of  the  effect  of  a  surface  pressure  traveling  with 
velocity  V  on  the  surface  of  an  elastic  half-space  (Reference  I3)  has 
been  one  of  the  most  useful  results  of  the  theoi^  of  elastic  wave  pro¬ 
pagation.  This  is  due  to  the  fact  that  this  problem  in  important  re¬ 
gions  represents  good  approximations  for  the  effects  of  air  blast  from 
an  atomic  weapon  in  the  vicinity  of  the  shock  front  (Figure  10).  The 
solution  (13)  has  been  a  major  tool  in  estimating  pressure  levels  cuad 
obtaining  shock  spectra  in  locations  at  some  distance  from  ground  zero. 

Obviously,  the  assumption  of  elastic  behavior  used  in  (I3)  is  an 
oversimplification  in  case  of  soil,  and  the  desirability  of  extending 
the  theory  to  allow  for  the  absorption  of  energy  is  quite  apparent.  The 
present  studies  are  intended  as  an  extension  of  (13)>  allowing  for  mech- 
anisiQs  of  energy  dissipation  to  be  expected  in  granular  materials. 

Under  the  current  contract,  as  reported  in  AFSWC-TIE-63-2  (Refer¬ 
ence  12)  a  megor  research  effort  was  made  to  introduce  two  types  of  dis¬ 
sipating  mechanisms  into  the  analysis,  generalizing  the  problem  treated 
in  (13)1  in  order  to  determine  the  effects  due  to  a  progressing  load  on 
the  surface  of  a  half -space  for  a  dissipating  material. 

One  type  of  material  considered  is  controlled  by  internal  friction 
of  the  Coulomb  type.  It  is  assumed  that  the  material  will  behave  line¬ 
arly  elastic,  undergoing  only  small  strains  if  the  internal  friction  is 
not  overcome.  If,  however,  certain  limiting  states  of  stress  etre  reached, 
internal  slip  will  occur.  Constitutive  equations  for  a  material  accord¬ 
ing  to  this  description,  i.e.  equations  which  give  the  relations  between 


stresses  and  strains  or  strain  rates,  can  te  formulated  on  the  basis  of 
plastic  potentials  In  the  manner  generally  used  in  the  analysis  of  elas¬ 
tic-plastic  bodies.  However,  when  forming  the  constitutive  equations  it 
vas  found  that  two  types  of  formulations  are  possible.  One,  proposed  in 
Reference  l4,  for  static  problems  applies  to  a  dilating  material,  l.e., 
to  a  material  which  Increases  In  volume  when  internal  slip  occiirs.  The 
other  formulation,  previoiuily  suggested  in  (13),  assumes  no  dilation  due 
to  slip.  Both  formilations  ere  considered  in  our  technical  report  (12). 

In  connection  with  the  formulation  of  the  boiixidary  condition  on  a 
free  surface  of  a  granular  medium,  it  vas  found  necessary  to  consider 
the  possibility  of  disintegration  of  the  medium.  The  situation  vas  ex¬ 
plored  on  a  mechanical  model,  simulating  a  medium  with  interior  friction 
of  the  Covilomb  type  (Appendix  B  of  12).  It  vas  concluded  that  such  dis¬ 
integration,  which  has  not  been  considered  previously,  occurs  when  a 
load  moves  with  subselsmlc  velocity  on  a  half-space  of  a  granular  medium. 
Such  disintegration  leads  to  dust  precursors  ahead  of  the  air  shock,  a 
matter  which  has  Indeed  been  observed. 

The  second  type  of  material  considered  in  which  energy  is  dis¬ 
sipated  is  locking  material.  Work  on  one  dimensional  wave  propaga¬ 
tion,  and  for  the  case  of  spherical  symmetry  (which  is  in  a  sense  also 
one  dimensional),  for  such  materials  was  discussed  earlier  in  this  report. 
The  extension  of  the  concept  of  locking  behavior  from  one  dlmenalon  to 
the  case  of  two  or  three  dimensions  is  a  major  step  of  a  fundamental 
nature,  and  the  result  obtained  under  the  contract  is  believed  to  be  a 
elgnlficant  achievement.  The  result  concerns  the  state  of  streas  and 
strain  at  the  front  of  consolidation.  In  the  one  dimensional  case  of 


l6 


propagation  of  plane  waves,  the  Rankine-Hugoniot  equations  give  a.  re¬ 
lation  between  the  stress  and  velocity  of  the  consolidation  front.  This 
simple  relation  suffices  to  formulate  the  appropriate  boundary  valvie 
problem.  In  two  or  three  dimensional  cases,  the  situation  is  much  more 
complicated,  as  the  mechanism  of  locking  in  three  dimensions  must  be 
introduced  before  the  conditions  at  the  locking  front  can  be  obtained. 

Once  the  conditions  on  the  consolidation  front  are  understood,  it 
is  possible  to  treat  two  dimensional  problems  for  locking  materials 
which  are  elastic  after  locking  or  have  other  behavior.  It  is,  there¬ 
fore,  possible  to  treat  the  locking  material  subject  to  Coulomb-slip 
after  consolidation,*  l.e.,  the  ceise  treated  in  (12)  can  be  coniDined 
with  locking  behavior. 

1,  Results  Obtained  for  Granular  MateriaLLs  Subject  to  Internal 
Coulomb  Friction 

For  this  material,  the  case  of  a  step  pressure  p  progressing  with 
a  velocity  V  on  the  surface  of  a  half-space.  Figure  11,  was  considered. 
Depending  on  the  state  of  stress,  the  material  described  above  will  be¬ 
have  elastically  without  slip,  or  its  strains  will  be  a  ccmobinatlon  of 
slip  and  elastic  effects,  such  that  the  solution  will  be  a  combination 
of  regions  with  slip,  and  without  slip.  To  simplify  the  approach,  and 
as  a  first  step  towards  a  more  complete  solution,  the  assiuaptlon  was 
made  that  in  regions  where  slip  occurs,  slip  strains  are  appreciably 
larger  than  elastic  strains,  such  that  the  latter  can  be  neglected. 
Based  on  this  assumption,  closed  form  solutions  for  the  problem  were 
obtained,  but,  because  of  the  above  eissumption,  they  are  restricted  to 


a  Bubseisiuic  range  of  the  velocity  V, 


V  ^  0.2  Cp  (1) 

vhere  c  ^  — _  is  the  velocity  of  p-waves  in  the  material, 

wuexB  Up  p(n.v)(i-2v) 

having  the  properties  p,  E,  and  of  the  non-slip  regions. 

The  character  of  the  solution  obtained  can  be  described  in  terms 
of  a  system  of  polar  coordinates  r,  the  origin  0  of  which  coincides 
with  the  front  of  the  applied  pressure  p.  Figure  11,  As  indicated  in 
this  figure,  there  are  three  wedge  shaped  regions,  separated  by  the 
lines  i  ■  45®,  i  «  135°.  Slip  occurs  only  in  one  of  them  45°  <  4  <  135® » 
in  the  two  others  the  material  acta  elastically.  The  expressions  for 
the  stresaes  are  given  in  (12).  It  is  important  that  the  vertical  stress 
for  i  o  0  (ahead  of  the  pressure  front  p)  does  in  general  not  vanish, 
but  that  in  this  location  particles  are  expelled  upwards  due  to  disinte¬ 
gration  of  the  material.  The  reactive  pressure  due  to  the  expelled  par¬ 
ticles  is  the  cause  of  the  vertical  stress  Just  below  the  surface  at 
i  a  0.  As  explained  in  detail  in  (iZ),  surface  disintegration  occurs 
for  all  values  of  the  applied  load  p  if  the  material  has  no  cohesion, 
but  only  for  -values  p  >  if  cohesion  is  present,  where  Pj^  is  a  limit¬ 
ing  -value  dependltig  on  the  value  of  cohesion. 

The  ultimate  purpose  of  the  analysis  being  the  design  of  under- 
groirnd  structures,  the  results  of  the  aimlysls  are  best  discussed  on 
two  pertinent  effects.  One  quantity  of  interest  is  the  total  permanent 
strain  at  a  point  due  to  the  passing  of  the  shock  wave  on  the  surface. 
This  strain  gives  an  indication  of  the  deformation  to  which  a  target. 
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say  a  cylindrical  shell,  Figure  12,  would  be  subjected  if,  in  first  ap¬ 
proximation,  the  effect  of  the  presence  of  the  target  on  the  free  field 
is  ignored.  For  a  typical  case,  k  =  I/3,  the  two  i>ermanent  principal 
strains  became 


c  =  o.in  (2) 

pr 

It  is  shown  in  (l,  p.  35)  that  in  the  range  of  velocities  given 
by  £q.  (l),  the  strains  (2)  are  appreciably  larger  than  in  an  elastic 
material  of  similar  properties  where  slip  does  not  occvir. 

A  second  quantity  of  interest  Is  the  acceleration  to  which  a 
target  Is  subjected,  llie  analysis  indicates  that  the  free  field  ac¬ 
celerations  *a"  are  at  each  instant  in  a  radial  direction  from  point  0  at 
the  head  of  the  applied  wave  towards  the  target.  The  value  "d' vanishes 
except  while  the  slip  region  passes  over  the  target.  The  tine  history 
of  the  acceleration  is  given  In  (l.  Figure  20).  Its  mnxinnm  value, 
for  k  =  1/3,  is 


max  a  =  0.24  -2.  (3) 

py 

where  y  is  the  depth  of  target.  Again,  it  is  demonstrated  In  (12)  that 
the  accelerations  In  the  material  with  slip  are  an  order  of  magnitude 
larger  than  in  an  elastic  material  of  similar  properties,  but  without 
slip. 

It  must  be  esq)hasized  that  the  results  obtained  In  (12)  are  not 
complete  and  further  work  is  required.  In  order  to  obtain  results  valid 


for  values  of  the  velocity  V  above  the  limit  in  Eq.  (l),  particularly 
in  -Che  superseismic  range,  it  is  necessary  to  extend  the  analysis  and 
include  the  presently  neglected  elastic  effects  in  the  slip  range.  This 
important  extension  is  planned  for  a  future  contract.  A  second  generali¬ 
zation,  involving  major  mathematical  difficulties,  concerns  the  treat¬ 
ment  of  decaying  applied  pressures,  instead  of  step  pressures.  In  the 
elastic  case,  the  differential  equations  are  linear,  such  that  the  treat¬ 
ment  of  decaying  pulses  is  easily  acconplished  by  superposition  of  the 
solutions  (13)  for  a  line  load.  In  the  present  case,  the  differential 
equations  are  non-linear,  such  that  the  decaying  pressure  wave  requires 
a  novel  analysis.  An  approach  to  this  problem  has  been  conceived,  and 
a  solution  will  be  attempted  in  the  future  contract. 

2.  Results  Obtained  Concerning  Twx?  Dimensional  Waves  in  Locking  Media 

The  mathematical  conditions  which  must  be  satisfied  on  the  locking 
front  for  the  case  of  one  dimensional  wave  propagation.  Figure  13,  are 
known  for  a  variety  of  types  of  material  behavior  after  locking,  i.e. 
rigid  or  linear  eleistlc.  These  conditions  are  derived  from  momentum  con¬ 
siderations  and  from  the  appropriate  one  dimensional  stress -strain  Iw 
in  a  relatively  simple  manner,  because  the  principal  stress  and  the  par¬ 
ticle  velocities  are  in  the  one  dlmei^sional  case  necessarily  in  the  di¬ 
rection  of  wave  propagation,  i.e.  at  right  angles  to  the  locking  front. 

The  situation  in  case  of  two  dimensional  wave  propagation.  Figure 
14,  is  much  more  cooiplex.  Considering  a  small  element  ds  of  the  consoli¬ 
dation  front.  Figure  15,  there  may  be  a  normal  stress,  0,  and  a  tangen¬ 
tial  one,  T .  Correspondingly,  the  velocity  of  the  particles  at  the  front 


may  have  noi’mal  and  tangential  components.  Using  roomentuia  considerations 
and  the  appropriate  stress-strain  law  for  plane  strain,  relations  between 
stresses,  particle  velocities^  and  the  velocity  of  propagation  of  the 
locking  front  can  be  formulated.  However,  the  formulation  becomes  coim- 
plicated  because  a  locking  material  is  a  non-linear  material,  and  the 
mathematical  description  of  general  non-linear  relations  between  stress 
and  strain  in  two  and  three  dimensions  requires  sophisticated  tools. 

In  spite  of  the  complicated  fomiulation  required  to  state  the 
problem,  quite  simple  relations  were  finally  obtained  for  isotropic 
locking  materials.  It  was  found  that  on  a  consolidation  front,  the  tan¬ 
gential  components  t  of  the  stress  must  vanish,  such  that  the  normal 
stress  0,  Flgui'e  15,  is  the  major  principal  stress.  As  a  consequence, 
the  particle  velocity  u  is  also  normal  to  the  front,  quite  similar  to  the 
case  of  one  dimensional  wave  propagation.  For  example,  when  the  lock¬ 
ing  strain  e  is  large  coiopared  to  subsequent  elastic  or  slip  strains, 
one  finds  a  relation  which  looks  jiist  as  the  result  for  one  dimensional 
locking  materials, 

■  w  M 

where,  however,  o  is  the  compressive  stress  normal  to  the  wave  front, 
end  U  la  the  normal  component  of  the  velocity  of  the  consolidation  front. 

In  addition  to  finding  that  the  shear  stresses  T  on  the  locking 
front  vanish,  such  that  the  normal  stress  0  is  a  principal  stress,  the 
analysis  also  furnishes  information  on  the  texigentlal  stresses  and  strains. 
It  is  found  tliat  the  tangential  strains  vanish.  In  the  case  of  a  locklxig 
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material  which  is  linearly  elastic  after  consolidation,  the  two  princi¬ 
pal  stresses  parallel  to  the  locking  front  are  therefore 


Alternatively,  if  the  material  after  locking  is  elastic,  hut  sub¬ 
ject  to  slip,  the  two  principal  stresses  parallel  to  the  locking  front 
are 

02  0  0^  =  ko  (5b) 

where  the  definition  of  k  (and  the  behavior  of  the  materieLL  in  general) 
are  as  previously  discussed  in  (12). 

Once  a  relation  between  a  and  U  of  the  type  of  Eq,  (4)  is  avail¬ 
able,  this  relation,  together  with  an  equation  of  the  type  (5)  for  the 
principal  stresses  Og  and  o^,  permits  the  formulation  of  a  boundary 
value  problem  for  the  region  behind  the  locking  front.  As  an  example, 
the  case  of  a  step  load  of  Intensity  p  progressing  with  velocity  V  on 
the  surface  of  a  half -space  has  been  tentatively  considered  for  the  case 
of  a  material  which  is  linearly  elastic  after  consolidation.  It  is  found 
that  there  is  a  plane  consolidation  front  making  an  angle  3  with  the 
swface.  Figure  l6,  where  p  dei)ends  on  the  values  of  p  and  V.  !Iho 
stresses,  velocities,  etc.  in  the  wedge  shaped  region  can  bo  obtained 
from  an  elastic  bouiidary  value  problem  for  this  region.  The  known  con¬ 
ditions  on  the  horizontal  surface  and  the  relatioiui  (4)  and  (5a)  on  the 
consolidation  front  furnish  the  proper  nutdber  of  boundary  conditions. 


Qlie  details  of  these  solutions,  particularly  for  the  case  of  laaterial 
BUhject  to  CouloBib  slip,  vill  be  considered  in  a  future  contract. 

Conciusiotis 

The  researches  summarized  in  this  report  show  that  air  induced  ground 
shock  effects  in  the  high  superseismic  range  can  be  approximated  by  means 
of  one  dimensional  wave  propagation  theory  in  a  non-linear  material  of 
decreasing  compressibility. 

Predictions  which  are  based  on  such  one  dimensional  calculations  show 
good  agreement  with  ground  motion  measurements  and  also  with  various 
currently  used  semi-empirical  methods. 

Numerous  additional  problems  which  arise  in  ground  motion  predictions 
and  in  the  evaluation  of  laboratory  experiments  require  the  solution  of 
two  dimensional  phenomena.  Under  this  contract,  results  were  also  ob¬ 
tained  for  a  step  pressure  progressiV\g  with  low  subseismic  velocity  over 
the  surface  of  a  half  space.  Two  distinct  models  are  considered.  In  one 
case  the  material  (depending  on  the  state  of  stress)  will  respond  elastically, 
but  at  other  states  slip  will  occur,  subject  to  internal  Coulomb  friction. 

In  the  second  case,  the  material,  exhibits  locking  characteristics  and  the 
propagation  of  a  locking  front  is  considered. 

The  results  of  these  researches  have  led  to  the  theoretical  proof  of 
the  existence  of  dust  precursors  which  have  previously  been  observed  but 
not  explained  in  these  terms.  Future  work  on  these  topics  should  consider 
the  use  of  a  pressure  pulse  advancing  with  superseismic  velocity  over  the 
surface  of  the  half  space. 


23 


REFERENCES 


1.  P.  Weidlinger,  "Analytical  and  Experimental  Studies  on  Locking  Media", 

29th  Syraposium  on  Shock,  Vibration  and  Associated  En¬ 
vironment,  Nov.  i960. 

2.  W.  Heierli,  "Die  Dynamik  eindimensionaler  Bodenkblrper  im  nichtlinearen, 

nichteleistischen  Bereich",  Doctoral  Thesis  E.T.H.  Zurich, 
Switzerland,  I96I. 

3.  A.  A.  Thon5)6on,  "A  Coo^jarison  of  the  Dynamic  and  Static  Stress-Strain 

(hirve  in  Sand  Ifader  Confined  and  Unconfined  Conditions", 

BRL  Memorandum  Report,  No.  1262,  April  I96O. 

k.  N.  M.  Newmark,  "Final  Report  of  Special  Studies  Related  to  the 
Wsio7A-2  Launcher  Installation",  AMF,  Jvine  20,  1958- 

5.  M.  Salvadori,  R.  Skalak  and  P.  Weidlinger,  "Waves  and  Shocks  in 

Liking  and  Dissipative  Media",  Proceedings  ASCE,  Vol.  86, 
No.  EM2,  April  19^). 

6.  M.  Salvadori,  R.  Skalak  and  P.  Weidlinger,  "Spherical  Waves  in  a 

Plastic  Locking  Medium",  Proceedings  ASCE,  Vol.  87,  No. 

EMI,  Feb.  1961. 

7.  R.  Skalalt  and  P.  Weidlinger,  "Attenxiation  of  Stress  Waves  in  Bi- 

Unear  Materials",  AFSWC-TN-60-30,  Oct.  i960. 

0.  S.  Sutcliffe,  "Strong  Shcok  Formation  in  Bi-Unear  Hardening  Media", 
Proceedings  ASCE,  Vol.  89,  No.  EM2,  April  1963. 

9.  P.  Weidlinger  and  A.  Matthews,  "A  Method  for  Prediction  of  Ground 
Shock  Phencxnena  in  Soils",  APSWC-TDR-6I-66,  Aug.  I96I. 

10.  ASCE,  "Design  of  Structures  to  Resist  Nuclear  Weapons  Effects",  ASCE 
Manual  of  Ekigineering  Preu:tlce,  No.  4-2. 

U.  M.  Salvadori  and  P.  Weidlinger,  "Induced  Ground  Shock  in  Granular 
Media",  STL,  Bulletin  Research  Division,  March  1958. 

12.  H.  H.  Bleich  and  E.  Heer,  "Step  Load  Moving  with  Low  Subseismic 

Velocity  on  the  Surface  of  a  Half -Space  of  Chranular  Mei- 
terial",  AFSWC-TDR-G3-2,  April  I963. 

13.  J.  Cole  and  J.  Hath,  "Stresses  Produced  in  a  HeLIf- Plane  by  Moving 

Loads",  Journal  of  Applied  Mech.,  Vol.  25,  p.443,  Dec.  I958, 


14.  D.  C.  Drucker  and  W.  Prager,  "Soil  Mechaxiics  and  Plastic  finalysis 

or  Limit  Design",  Quart.  App!’ .  Math.,  P.  15'j',  I952. 

15.  A.  Y.  Ishlinski,  "On  Two  lUmensional  Motion  of  Sand"  (in  Russian), 

Ukranian  Math.  Journal,  Inst,  of  Math.,  Vbl.  VI,  Ko.  k, 

195^. 


VELOCITY  V 


F  I  G  .  I  I 


F  I  G. I  2 


» 


23 


VELOCITY  V 


CONSOLIDATION  FRONT 


VELOCITY  U 


F I G. 16 


29 


TDR-63-3048 


No,  Cvs 


2 

1 

1 

1 

1 


1 

n 

2 


1 

3 

3 

2 

1 

1 

2 

2 


1 

5 

30 

5 


DISTRIBUTION 

HEADQUARTERS  USAF 
Hq  USAF  ( AFOCE),  Wash,  DC  20330 
Hq  USAF  (AFRNE-B,  Maj  Lowry),  Wash,  DC  20330 
Hq  USAF  (  AFTAC),  Wash,  DC  20330 

USAF  Dep,  The  Inspector  General  ( AFTDI) ,  Norton  AFB,  Calif 
9  2409 

USAF  Directorate  of  Nuclear  Safety  ( AFINS) ,  Kirtland  AFB,  NM 
87117 

MAJOR  AIR  COMMANDS 
AFSC  (SCT),  Andrews  AFB,  Wash,  DC  20331 
AUL,  Maxwell  AFB,  Ala  36112 

USAFIT  (USAF  Institute  of  Technology),  Wright-P atter son 
AFB,  Ohio  45443 

AFSC  ORGANIZATIONS 

ASD  ( ASNXRR,  Tech  Info  Ref  Br.  ,  Reports  Div.  ,  Marie  Koeker), 
Wright-Patterson  AFB,  Ohio  45433 

BSD(BSR),  Norton  AFB,  Calif  92409 

SSD(SSSC),  AF  Unit  Post  Office,  Los  Angeles,  Calif  90045 
ESD(ESAT),  Hanscom  Fid,  Bedford,  Mass  01731 
AF  Msl  Dev  Cen  (RRRT),  Holloman  AFB,  NM  88330 
AFMTC  (MU-135,  ATTN:  MTBAT),  Patrick  AFB,  Fla  32925 
APGC(PGAPl),  Eglin  AFB,  Fla  32542 
AEDC  (AEOI),  Arnold  Air  Force  Station,  Tenn  37389 

KIRTLAND  AFB  ORGANIZATIONS 
AFSWC,  Kirtland  AFB,  NM  87117 
(SWEH) 

(SWT) 

AFWL,  Kirtland  AFB,  NM 
(  WLL) 

(WLRS) 


30 


TDR  6  3-3048 


No,  Cvs. 


2 

1 

2 

1 

1 

2 

1 


1 

1 

1 

1 

1 

1 


1 

2 


DISTRIBUTION  (cont'ri) 


OTHER  AIR  FORCE  AGENCIES 

Director,  USAF  Project  RAND,  via:  Air  Force  Liaison  Office , 
The  RAND  Corporation,  1700  Main  Street,  Santa  Monica,  Calif 
90406 

AFOAR,  Bldg  T-D,  Wash,  DC  20333 
ARMY  ACTIVITIES 

Director,  Ballistic  Research  Laboratorie  s  (  Library) , 

Aberdeen  Proving  Ground,  Md  21005 

Hq  US  Arnay  Air  Defense  Command  ( ADGCB) ,  Ent  AFB,  Colo 
80912 

Director,  Army  Research  Office,  3045  Columbia  Pike, 

Arlington,  Va  22204 

Director,  US  Army  Waterways  Experiment  Sta  (  WE  SRL), 

P.  O.  Box  631,  Vicksburg,  Miss  39181 

Director,  US  Army  Engineer  Research  &  Development 
Laboratories,  ATTN:  Tech  Documents  Center,  Ft  Belvoir,  Va 

NAVY  ACTIVITIES 

Chief,  Bureau  of  Yards  and  Docks,  Department  of  the  Navy, 

Code  22.  102,  Wash  25,  DC 

Commanding  Officer  and  Director,  David  Taylor  Model  Basin, 
Wash  7,  DC 

Commanding  Officer  and  Director,  Naval  Civil  Engineering 
Laboratory,  Port  Hueneme,  Calif 

Commander,  Naval  Ordnance  Test  Station  (  Code  12), 

China  Lake,  Calif  93557 

Office  of  Naval  Research,  Wash  25,  DC 

US  Naval  Weapons  Evaluation  Facility  {NWEF)(Code  404), 

Kirtland  AFB,  NM  87117 

OTHER  DOD  ACTIVITIES 

Chief,  Defense  Atomic  Support  Agency  { Document  Library,  Miss 
Gertrude  Camp),  Wash,  DC  20301 

Commander,  Field  Command,  Defense  Atomic  Support  Agency 
{FCAG3,  Special  Weapons  Publication  Distribution),  Sandia  Base, 
NM  87115 


31 


TDR  63-3048 


No.  Cvs. 
1 

1 

20 


1 

1 


1 

1 

2 

1 

2 

2 

2 

2 

1 

1 

2 

1 


DISTRIBUTION  (cont'd) 


Director,  Advanced  Research  Projects  Agency,  Department  of 
Defense,  The  Pentagon,  Wash,  DC  20301 

Office  of  Director  of  Defense  Research  &  Engineering,  ATTN: 
JohnE.  Jackson,  Office  of  Atomic  Programs,  Rm  3E  1071, 

The  Pentagon,  Wash  25,  DC 

Hq  Defense  Documentation  Center  for  Scientific  and  Technical 
Information  ( DDC) ,  Bldg  5,  Cameron  Sta,  Alexandria,  Va  22314 

AEC  ACTIVITIES 

Sandia  Corporation  { Information  Distribution  Division),  Box  5800, 
Sandia  Base,  NM  87115 

Sandia  Corporation  ( Technical  Library),  P.  O.  Box  969, 
Livermore,  Calif  94551 

OTHER 

OTS,  Department  of  Commerce,  Wash  25,  DC 

Space  Technology  Labs,  Inc.,  ATTN;  Information  Center, 
Document  Procurement,  P.  O.  Box  95001,  Los  Angeles  45, 

Calif 

University  of  Illinois,  Talbot  Laboratory,  Room  207, 

Urbana,  Ill 

Massachusetts  Institute  of  Technology,  Dept  of  Civil  &  Sanitary 
Engineering,  77  Massachusetts  Avenue,  Cambridge,  Mass.  , 
ATTN;  Prof.  R.  V.  Whitman 

South  Dakota  School  of  Mines  and  Technology  (ATTN;  Prof. 
Edwin  H.  Oshier),  Rapid  City,  SD 

AF  Shock  Tube  Facility,  Box  188,  University  Station, 
Albuquerque,  NM 

IIT  Research  Institute,  ATTN;  Dr.  Sevin,  3422  South  Dearborn 
Street,  Chicago  15,  Ill 

MDR  Division,  General  American  Transportation  Corp.  , 

7501  North  Natchez  Ave.  ,  Niles  48,  111 

National  Engineering  Science  Co.  ,  ATTN;  Dr,  Soldate, 

711  South  Fairoaks  Ave.  ,  Pasadena,  Calif 

United  Research  Services,  1811  Trousdale  Drive,  Burlingame, 
Calif 

Stanford  Research  Institute,  ATTN;  G-037,  External  Reports, 
Menlo  Park,  Calif  94025 

University  of  Notre  Dame,  Dept  of  Civil  Engineering,  ATTN: 

Dr.  Saxe,  Notre  Dame,  Ind 


32 


TDR  63-3048 


No.  Cvs. 
1 

1 

2 

1 

20 

1 

1 

1 

2 

1 

1 

1 

1 

1 

1 

1 

1 

1 


DISTRIBUTION  (cont'd) 


Stanford  University,  School  of  Mechanical  Engineering,  ATTN: 
Dr.  Jacobsen,  Stanford,  Calif 

University  of  Washington,  ATTN;  Dr.  I.  M.  Fyfe,  Seattle  5, 
Wash 

Purdue  University  School  of  Civil  Engineering,  ATTN:  Prof. 
Leonards,  Lafayette,  Ind 

California  Institute  of  Technology,  ATTN:  Dr.  Paul  Blatz, 

1201  East  California  Blvd,  Pasadena,  Calif 

Paul  Weidlinger  and  Associates,  770  Lexington  A.ve,  New  York 
21,  NY 

Shannon  and  Wilson,  Foundation  Engineers,  1105  North  35th 
Street,  Seattle  3,  Wash,  ATTN:  S.  D.  Wilson 

Mitre  Corporation,  Bedford,  Mass 

University  of  California,  Berkeley,  Calif.  ,  ATTN:  Prof. 

H.  B.  Seed 

University  of  Michigan,  Ann  Arbor,  Mich,  ATTN;  Prof. 

F.  E.  Richart  and  Prof.  Bruce  G.  Johnston 

RAND  Corporation,  1700  Main  Street,  Santa  Monica,  Calif 
90406  ATTN;  H.  L.  Brode 

Ralph  M.  Parsons  Co.  ,  6  17  West  7th  Street,  Los  Angeles  17, 
Calif,  ATTN:  Herbert  R.  Saffell 

Colorado  School  of  Mines,  Golden,  Colo.  ,  ATTN:  David  Card 

Physics  International,  2229  Fourth  Street,  Berkeley  10,  Calif 
ATTN:  Dr.  C.  S.  Godfrey 

Northrup- Ventura,  1515  Rancho  Conejo  Blvd,  Newbury  Park, 
Calif.,  ATTN;  Dr.  JohnG.  Trulio 

Structural  Behavior  Engineering  Laboratories,  3515  West 
Clarendon  Avenue,  Phoenix  19,  Ariz,  ATTN:  Mr.  J.  D.  Hess 

St.  Louis  University,  Institute  of  Technology,  ATTN;  Dr.  Carl 
Kisslinger,  36  21  Olive  Street,  St.  Louis  8,  Mo 

Northwestern  University,  Evanston,  111.  ,  ATTN:  Dr.  Robert 
L.  Kondner 

Official  Record  Copy(WLRS,  Mr.  Walsh) 


33 


